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Abstract 

The linear and nonlinear properties of the Rao-dust-magnetohydrodynamic (R-D-MHD) waves 
in a dusty magnetoplasma are studied. By employing the inertialess electron equation of mo- 
tion, inertial ion equation of motion, Ampere's law, Faraday's law, and the continuity equation 
in a plasma with immobile charged dust grains, the linear and nonlinear propagation of two- 
dimensional R-D-MHD waves are investigated. In the linear regime, the existence of immobile 
dust grains produces the Rao cutoff frequency, which is proportional to the dust charge density 
and the ion gyrofrequency. On the other hand, the dynamics of an amplitude modulated R-D-MHD 
waves is governed by the cubic nonlinear Schrodinger equation. The latter has been derived by 
using the reductive perturbation technique and the two-timescale analysis which accounts for the 
harmonic generation nonlinearity in plasmas. The stability of the modulated wave envelope against 
non-resonant perturbations is studied. Finally, the possibility of localized envelope excitations is 
discussed. 

PACS numbers: 52.27.Lw, 52.35.Hr, 52.35.Mw, 52.35. Sb 



* Preprint submitted to Physics of Plasmas. 



1 



I. INTRODUCTION 



A wide variety of electrostatic and electromagnetic oscillatory modes are known to prop- 
agate in unmagnetized and magnetized plasmas Since more than a decade ago, it 
has been pointed out, and is now well established, that the presence of heavy charged dust 
particulates in a plasma may strongly modify the dispersion properties of the known low- 
frequency modes, and may also introduce novel waves 3,0]. For instance, inclusion of the 
dust particle dynamics in an unmagnetized dusty plasma gives rise to the dust-acoustic 
waves , while the modification of the plasma constituents' charge balance is responsible 
for the dust ion-acoustic waves |^ , characterized by an increased phase speed in comparison 
with the acoustic speed in an electron-ion plasma without dust. In a magnetized dusty 
Dlasma, a variety of new modes have been shown to exist, including modified Alfven waves 
7| propagating alongto the direction of the external magnetic field B, as well as the modified 
magnetoacoustic fl and drift-electromagnetic Q waves propagating across B. 

In this paper, we will focus on the linear and nonlinear properties of the Rao-dust- 
magnetohydrodynamic (R-D-MHD) waves in two space dimensions. The dispersion 
characteristics of the two-dimensional (2D) R-D-MHD waves differ from the ordinary mag- 
netosonic waves propagating in a magnetized electron-ion (e-i) plasma; of particular im- 
portance is the existence of a novel cutoff frequency due to the presence of charged dust 
grains, as first reported by Rao in his classic paper 9]. Apart from being interesting from a 
fundamental point of view, and not so widely studied so far, the R-D-MHD waves have been 
recently shown fll^ to be excited by the upper-hybrid waves in a uniform dusty magneto- 
plasma. Our objective here is twofold: i) to present two-dimensional R-D-MHD modes, ii) 
to study the amplitude modulation of finite amplitude 2D R-D-MHD waves. Assuming the 
existence of a uniform external magnetic field and relying on the two-fluid model description, 
we will calculate analytically the harmonic response of the system to a small displacement 
from equilibrium, trying to point out the role of the dust. The nonlinear modulation of 
the wave's amplitude will then be considered by making use of an appropriate reductive 
perturbation method jl^ . Iisl . Q| . The R-D-MHD wave stability will then be investigated 
and the existence of envelope excitations will be discussed. 

The manuscript is organized in the following fashion. In Sec. II, we present the govern- 
ing equations for the R-D-MHD waves. Linearized equations and harmonic solutions are 
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presented in Sec. III. Considering oblique nonlinear amplitude modulations of finite ampli- 
tude R-D-MHD waves, we derive the nonfinear Schrodinger equation in Sec. IV. A stabifity 
analysis is carried out in Sec. V. Section VI contains a discussion of localized R-D-MHD 
modes. Our conclusions are highlighted in Sec. VII. 



II. THE MODEL 



We consider a three-component fully ionized dusty plasma composed of electrons (mass 
m, charge e), ions (mass mj, charge qi — +Zie) and heavy charged dust particulates (mass 
TTid, charge qa — sZ^e), henceforth denoted by e, i, d respectively. Dust mass and charge 
will be taken to be constant, for simplicity. Note that both negative and positive dust charge 

cases are considered, distinguished by the charge sign s = sgn qd = ±1. 

The plasma is immersed in a uniform external magnetic field along the 5— direction: 
Bo = -Bo i (-Bo = const.) 



A. Evolution equations 

Let us consider the MHD system of equations for electrons and ions. The massive dust 
particles are assumed to be practically immobile ('frozen' i.e. ~ n^ o), since we are inter- 
ested on timescales much shorter than the dust plasma period (~ The electron/ion 
number density rij^e and velocity Vj_e a-re governed by the equations 

drif, 



dt 
'dt 



+ V-(neUe) = 0, (1) 

+ V-(niU,) = 0, (2) 



E + X B = , (3) 

and 

rrii Di Ui = rrii + Ui • Vui^ = e (^E + -^u; x B^ , (4) 

where we have completely ignored the electron inertia, as well as pressure (temperature) 
effects (for all species a); the convective derivative operator: Di = ^ + Ui • V has been 
defined. E and B denote the (total) electric and magnetic fields, E = -|- Ei and B = Bq -|- 
Bi, respectively, i.e. index (1) denotes the external (wave) field components. Throughout 
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this text, we shall assume that Ei = {Ei^^, -E-i.y, 0) and Bi = (0, 0, Bi), where Ei^^/y and Bi 
are allowed to depend on {x, y, t). The system is closed with Maxwell's equations; neglecting 
the displacement current, Ampere's law reads 

VxB = — J = — Vgari„Uc,= [Zi rijUi - ngUe) (5) 

c c V c 

and Faraday's law is 

Note that the condition V ■ B = here reduces to dB jdz = 0. At equilibrium, the overall 
neutrality condition holds 

Uefi - Zi Tlifl - sZdUd = 0. (7) 

Since we are interested in waves propagating in the direction perpendicular to the mag- 
netic field, will shall assume, throughout this study, that the velocities Uq, (a = e,i), the 
wavenumber k and the electric field E lie in the xy— plane. See that E is orthogonal to Ue 
and B, due to (jHI- 



B. Reduced system of equations 



By eliminating E in (jHI), (jH), we obtain 



rrii Di Ui = Zi - (ui - Ue) xB 

c 



which, combined with (jSl), in order to eliminate Ue, i.e. 



(8) 



Up 



^U; - (V X B) 



yields 



mi Di Ui = Zi ^ (ui X B) + ( V X B) X B 

ripC Anrip 



Z.^(uiXB) 



B ■ VB - - VB^ 
2 



(9) 



(10) 



where we have used the quasineutrality condition, rig — ZiUi — sZdUd = 0. We observe 
that, to a first approximation, i.e. assuming very weak magnetic field non-uniformity, the 
ions (and the electrons due to 0) are subjected to a rotation due to the presence of charged 
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dust grains, as also shown in Ref. 10|: notice the Lorentz centripetal force in the right- 
hand-side of (fTUj) . associated with a rotation frequency which is directly proportional to the 
dust charge (and vanishes without it). 

Now, by eliminating E in (jH)), (jHI) and using Q, we obtain 



OB 
'dt 



Vx 



ZiTli 



Up 



Ui X B) 



47re 



Vx 



— (V X B)xB 



Hp 



(11) 



Note that Eqs. (jH)) - (fTT|) lead to a novel low-frequency electromagnetic mode, associated 
with the presence of charged dust grains, as was recently shown in Ref. cf. Eqs. (4), 
(6) - (8) therein. 

The system of equations (fTU]). (|TT|) is not closed in B and Ui, since it also involves Up and 
Hi (both variable), unless one limits the analysis to small (first order) perturbations from 
equilibrium. Otherwise, for a consistent description, one should either use the complete 
system of Eqs. (HI) - (jH)) or retain Eqs. ((H), 0, (0), (0) and (jHI) instead. In the following, 
we will adopt the former option. 

The set of equations (H)) to (0) is a closed system describing the evolution of the state 
vector S = (ng, n^, Ue, Ui, E, B). By assuming that no other vector quantity has a compo- 
nent along the magnetic field B = B z = {Bo + Bi) z, viz. E = + Ei = [E^, Ey, 0), and 
Ue/i = {ue/i,x, '^e/i,y, 0), whcrc E^/y, u^/y and Bi are functions of {x, y; t}, we obtain 



+ — {Ue Up^^) + — {Up Up^y) = , 

dt dx dy 

dt dx dy 
1 



rrti 



d_ 
di 
d 



d_ 

OX 



d_ 

dy 



Ui 



Ey = +-Up^^B , 
c 

= ZiC (e^ + -Ui^yB 



d 



d 



dt 



' dx 



dy 



mi — + Ui^x-^ + Ui,y-^ )Ui^y = ZiCiEy Ui^^B 



ZicB 
c 

ZicB 
c 



{Ui^y Up^y) , 



[Ui 



dB _ 47re 

dy 
dB _ 
dx c 



c 

Ane 



(^ZifliUi^x ^e^e,a;) ; 
ZiTLiUi^y ^e^e,y) 5 



(12) 

(13) 
(14) 
(15) 
(16) 
(17) 
(18) 
(19) 
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and 

dx dy c dt ' 
describing the evolution of the 9 scalar quantities: n^, rii, Ue^x/y, Ui^^/y, E^/y and B. Note 
that ()14|) . ()15p can be used to eliminate E in (PUI) . which then becomes 

dB _ d{u,,xB) d{ue,yB) 



dt dx dy 

Equations ()12|) - ()20|) will be the basis of the analysis that follows. 



(21) 



III. LINEARIZED EQUATIONS - HARMONIC SOLUTIONS 

By linearizing around the equilibrium state So = (?^e,o, ^i.o, 0, 0, 0, Bq) viz. S = So + Si 
and assuming linear perturbations of the form: Si = Si exp z(kx— tut) +c.c. = Si expi(/ca;+ 
ky — ut) + c.c. ('c.c' denotes the complex conjuguate) we obtain a new system of (linear) 
equations for the perturbation amplitudes (5'i)^.. A tedious, yet perfectly straightforward 
(see in the Appendix), calculation leads to the equations 

u{iuvx + 5Vtc,iVy) = iQl^L'^kxiKvx + kyVy) , 

u{iuvy - SQc,iVx) = iVtliL'^ky{kxVx + kyVy) , (22) 

in terms of the ion velocity component amplitudes Vj = unj {j = x,y), where we have 
defined 

- the ion gyrofrequency: fi^j = , 

- the characteristic length: L = ( 4^'eC^'° ) ' 

- the (dimensionless) dust parameter: S = -^-^ = s ^1 — -^-^^ ; see that 5 cancels in 
the dust-free limit [cf. d?))]. 

Equations (P^ . b) constitute a 2 x 2 homogeneous Cramer (linear) system, in terms of 
Uy, whose determinant should vanish in order for a non-trivial solution to exist; the wave 
frequency Q and wavenumber k are thus found to obey the dispersion relation 



UJ 



+ C^e (23) 



g 

where k = {k^ + k^Y^'^] we have defined 
- the 'gap frequency' Ug 



Zd^dfiZieBQ , . 

ujg = ■ = 6 i2c,i (24) 

riefirriiC 
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and 

- the characteristic velocity C = Qc,i L, given by 



^2 _ ZfE^rLifl 



ZiU. 



i'H,0 



Bl 



'l-s5fVl 



(25) 



i.e. C = ^c'iL = (1 — s5) Va, where Va = Bq/ {Aixni ^raiY^'^ is the Alfven speed. Notice the 
effect of the dust, which results in 

- a finite {'gap') oscillation frequency at the infinite wavelength {k —>■ 0) limit, and 

- a modified phase speed Vph = oj/k Vg = C^k/u, for 6 ^ 0); as a matter of fact, the 
phase speed Vg (~ C for lj ^ ujg) is higher (lower) than the Alfven speed Va in the presence 
of negative (positive) dust. 

n 

Notice that (j2Sl) coincides with (9) in Ref. [10|. It should also be pointed out that the 
existence of both the cutoff frequency ujg and the modified Alfven speed C, associated with 
;he dust-magnetosonic waves, was predicted for the first time by Rao in his classic paper 



a 



The harmonic perturbation amplitudes Sij may now be calculated. Assuming k 
{kcos6, ksin6), one obtains the following relations 



and 



Hi 



^efi A _ (11) A 

— — i^i = Ci i^i , 

-Do 



B, ^ cf) B, 



to COS 9 — i^l^j 



ZiBo 



U 



uj sin 6 + iQ 



-1 



ZjU. 



sin 6 



cos( 



^1 _ Jii) h 

= C3 -Dl 



Bok 
Bi 

Bnk 



Ui 



^^''^ [u COS 9 + i sSflc i sin 9] — \- = Cg^^'' Bi 
ZiUifl ' Bok 



u. 



[u sm 9—1 sd\lc,i cos 



Bn 



-U 



Ce Bi , 



' Bok 
cos 6* + sin 6* 



Bi 
ck 



c^-r'^B^ (32) 



(26) 
(27) 
(28) 
(29) 
(30) 
(31) 



c 



iQ 



-1 

c,i 



Bi 



s\Yi9 - UJ cos9\^ = cl"^ Bi (33) 
ck 
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(obviously, Cg =1). Note that these relations satisfy 

El ■ Ue,l = , 



in agreement with Q; also, 

k ■ Ue,i = k ■ Ui 1 = —Bi 

ne,0 -DO 

as expected (see the Appendix) as well as 

(remember that the amplitudes Ue/i^x/y are complex) implying that the ions and electrons 
oscillate in (out of) phase for uo lower (higher) than $7^,4; i-e. for wavenumber values k below 
(above) a threshold k^r = ^^^^^-^(Y^fl)^^^ (see that k^r — > in the case of complete electron 
depletion in the plasma, i.e. 5 — > 1, s = —1). 



IV. OBLIQUE NONLINEAR AMPLITUDE MODULATION 

Let us consider the system (jl2j) - (j2(Jj) . which describes the evolution of the (nine scalar) 
components of S: {n^, rii, Ue,x, Ue,y] Ui^x, Ui^y] Ex, Ey] Bi}. 

In order to study the amplitude modulation of the R-D-MHD waves presented in 
the previous section, we will assume small deviations from the equilibrium state S*^*^-* = 
{ne,o, rii^o] 0, 0; 0, 0; 0, 0; Bq] by taking 

oo 

S = + e S(i) + S(2) + ... = + E S^") , 

n=l 

where e -C 1 is a smallness par ameter. Following the standard multiple scale (reductive 
perturbation) technique [l^ Il3|. we shall consider the stretched (slow) space and time 
variables 

( = e(a; - At) , r = e^t, (34) 

where A, having dimensions of velocity, is a real parameter to be later defined. In order to 
allow for an oblique amplitude modulation on the R-D-MHD wave, we will assume that all 
perturbed states depend on the fast scales via the phase 6i = k • r — cut = kxX + kyy — ut 
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only, while the slow scales enter the argument of the /— th harmonic amplitude Si"'\ allowed 
to vary only along x, 

oo 
l=—oo 

The reality condition S^/ = sl""* is met by all state variables. Note that the (choice of) 
direction of the propagation remains arbitrary, yet modulation is allowed to take place in 
an oblique direction, characterized by the angle variable 9. Accordingly, the wave-number 
vector k is taken to be k = (fc^:, ky) = {k cos 6', k sin 6*). According to these considerations, 
the derivative operators in the above equations are treated as follows 

d d d 2 ^ 



and 



■^-[V. + e—, 



V + ex— = (y, + e—,Vy), 



i.e. explicitly 

V.A'r'^e'''' ={Uk cos 9 A'i^^ +ex^^)e'''\ 

and 

V,aJ") e^''^ =Uk sin 9 A^""^ e'''\ 

for any of the components A^j^^ {j = 1, 9) of s["''. 

By substituting the above expressions into Eqs. (fT^ - (pUjl and isolating distinct orders in 
e, we obtain a set of (nine) reduced equations at each (nth-) order, describing the evolution of 
the (nine) components of S^"-*. The system is then solved (for each harmonic /), substituted 
into the subsequent order, and so forth. This is a rather standard procedure in the reductive 
perturbation method framework Q], and we shall not burden the presentation 

with unnecessary details. The outcome of the long algebraic calculation is presented in the 
following, while essential details are presented in the Appendix. 

The first order [n = 1) first harmonic (/ = 1) equations are just as described in the 
previous section. Recall the (parabolic) form of the dispersion relation (jSHI), which arises as 
a compatibility condition. The amplitudes of the first harmonics of the perturbation, say 
^ij (j = •••jQ) (i.e. precisely Aj^i in the previous Section), then come out to be directly 
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proportional to the magnetic field perturbation, viz. A^^j = c^j^^^b[^^; the coefficients c^-^^^ 
are defined in - ()33|) above. Only the first harmonics have a contribution at this order; 
indeed, for n = 1,1 = 0, one obtains a (6 x 6) linear homogeneous system of equations for the 
(6 components of) Ue, Ui, E; interestingly, the determinant D^^ ~ ^c,i(ld is non-zero due to 
(and only in) the presence of dust, so we obtain the trivial solution for the zeroth-harmonic 
contribution, u^^q = u^- = Eq^-* = 0. In addition, n^^l = = Sq^-* = 0, as imposed by the 
(n = 2, Z = 0) equations. 

A. Second order in e: group velocity, 0th and 2nd harmonics 

The second order {n = 2) equations for the first harmonics provide the compatibility con- 
dition: A = duj /dkx = uj'ik) cos 6*, which defines A as the group velocity Vg = {C^k/u) cos9 
(the characteristic velocity C was defined previously). The 2nd-order corrections to the first 
harmonic amplitudes, say A^j {j = 1, 9), come out to be A^j = cf^^ dB[^^ /d(, where the 
coefficients c^^^^ are presented in the Appendix. 

As expected, second order harmonic contributions arise in this order; their amplitudes, 
defined by the equations for n = 2, / = 2, are found to be proportional to the square of the 
first order elements, e.g. in terms of b[^^: A'^j = c^^'^\B[^^y . The nonlinear self-interaction 
of the carrier wave also results in the creation of a zeroth harmonic, to this order; its strength 
is analytically determined by taking into account the / = component of the 3rd and 4th 
order reduced equations. The result is conveniently expressed in terms of the square modulus 
of the {n = 1,1 = 1) quantities, e.g. in terms of \b[^^\^ = (5^ ^)* b[^\ viz. = cf \b[^^\^ 
{j = 1, 9); once more, the definitions of c^-^^\ c^-^^"* can be found in the Appendix. Notice 
(see the Appendix) the dependence of the expressions derived in this Section (except those 
for He^i, B, in fact) on the value of 6. 

B. Derivation of the Nonlinear Schrodinger Equation 

Proceeding to the third order in e (n = 3), the equation for Z = 1 yields an explicit 
compatibility condition to be imposed in the right-hand side of the evolution equations 
which, given the expressions derived previously, can be cast into the form of the Nonlinear 
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Schrodinger Equation (NLSE) 

^|^ + ^^ + ^?HV = 0. (35) 

where i/j = b[^^ denotes the amphtude of the first-order electric field perturbation. Recall 
that the 'slow' variables {C,t} were defined in (j34j) . 

The dispersion coefficient P is related to the curvature of the dispersion curve as P = 

; the exact form of P reads 



1 ^ 

2 dkl 



sin^ d 



LU" (k) COS^e +Uj'{k) ^ 



P{k) = ^ (^cjj cos^ 9 + sin^ 9j , (36) 

which is positive for all values of the angle 6, as expected from the parabolic form of u;(k). 
The nonlinearity coefficient Q is due to the carrier wave self-interaction. It is given by 



Q 



-326^^ 77^ Zi {uefi - ZiUiflY {nefi + ZiUifi) + (? k^ mi TT (ne,o - ZiUi^ofiri^^Q - Zln\Q 



(37) 



Quite surprisingly, Q comes out to be independent of the angle d. However, as expected, 
the presence of charged dust grains in the charge balance equation strongly affects the 
numerical value of Q] notice, in passing, that this expression is not valid in the absence of 
dust grains (since the denominator then vanishes). 

The last expression for Q can be conveniently re-arranged, by making use of appropriate 
plasma quantities. Let us first define the dust parameter: = Uefl/ {ZiUifi)] see that: 
11 = 1 + s (ZdUdfi)/ (ZiUifi), due to ((Zj), so a value lower/higher than 1 corresponds to 
negative/positive dust charge sign; fi obviously tends to unity in the absence of dust (in any 
case, /i > 0). Check that /i = (1 — s6)~^ [or 6 = s(l — l/fi)], where 6 was defined above. 
By normalizing the wavenumber k as k = KtUp^i/c = {AixUiflZie^ /miC^Y^'^K {cOp^i is the ion 
plasma frequency), expression (|H7j) can be cast into an elegant form 

[{fi-if + x^Y^^zy 



(38) 



fi{fi-iymjc^ Q,^i 

(flc,i denotes the ion gyrofrequency defined previously). Retaining the approximate long- 
wavelength (i.e. vanishing wavenumber) behaviour of Q, we have 

g^K«l,,).-2^1^4^^ (39) 
11 



which is always negative and thus ensures, as we shall see in the following, stability at long 
wavelengths. Note, for later reference, that the same scaling results in relations ()23|) and 
(IHU)) taking, respectively, the reduced forms 



00 



and 



P 



X 



2-, 1/2 



(40) 



+ — sin^ e 



2 -1 3/2 



(41) 



V. STABILITY ANALYSIS 



The modulational stability profile of a carrier wave whose amplitude is described by the 



NLS 
here 



5quation 


m 




las 


m 


14, 


15 


iq 


121 



las long been studied, so only the main results have to be summarized 



The analysis consists in considering the linear stability of the monochromatic (Stokes's 
wave) solution of the NLSE ^ ip = {jj e'^\^\^^ + c.c. If the product PQ of the NLS 
coefficients is positive, the wave's envelope may develop an instability when subject to an 
external perturbation characterized by a wavenumber k lower than fccr = \J'^^\i^o\- The 
instability growth rate cr = \Imuj{k)\ then reaches its maximum value for k = kcr/V^, viz. 
<^max = IQI iV'oP- On the other hand, the wave will be stable for all values of k if the product 
PQ is negative. 

In our case, the dispersion coefficient P is positive, so one need only investigate the sign 
of the nonlinearity coefficient Q, which is entirely determined by the quantity in brackets in 
the right-hand-side of (jHEl); this is in fact a bi-quadratic polynomial of K, say p{K,fi). It is 
a matter of straightforward algebra to show (and an easy matter to confirm, numerically) 
that p{K,fi) (and Q) is negative for values of fi below ficr = 25.1146, i.e. for all values of 
the wavenumber x. Therefore, for negative dust charge [s = —1 i.e. /i < 1), the wave will 
always be stable. On the other hand, for positive dust charge (s = +1 i.e. /i > 1), the wave 
may become unstable (only) for values of /i above ficr i-e. in the case of positive dust charge 
concentration qdrid higher than ^ 2Aqini (a very rare situation, physically speaking, which 
implies a very high ion depletion in the plasma). The numerical value of Q, as expressed 
by relation (jH^ . is roughly depicted in figure [T] versus the wavenumber K and the dust 
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parameter /i. As predicted above, Q (and PQ) only reaches positive values for /i beyond 
^ 25 and K above ~ 10 (i.e. k > lOup^i/c), which is a hardly ever realizable physical 
situation. We conclude that the R-D-MHD waves are modulationally stable, in the presence 
of negatively charged dust grains, and (practically) also for positively charged ones. 



VI. LOCALIZED MODES 



Different types of envelope excitations (solitons) are known to satisfy Eq. ()H5|1 : in specific, 



one finds bright- (dark- or grey-) type so^ 



value of the coefficient product PQ |l3 . 



optics 



itons, e.g. pulses (holes) for a positive (negative) 
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16l Il7j , as already long known from nonlinear 



19|. According to the conclusions of the preceding Section, the R-D-MHD waves 



considered in this study will (in the majority of physically realizable situations) rather favour 
dark-type localized excitations, i.e. field dips (voids) propagating at a constant profile, 
thanks to the balance between the wave dispersion and nonlinearity. The analytical form of 
these excitations, depicted in Fig. |21 reads ijj{(,T) = J p{(,t) e^^^'^''^\ where 



P = Po 



a sech^ 



c 



U T 



(42) 



where a is a real parameter measuring the depth of the field void: < a < 1 {a = 1) 
corresponds to grey (black) solitons; see Fig. (|2b)- The complex expressions for the 
parameters a and G in the above expression fas well as related ones for bright solitons) can 



readily be found in the references 



xpression fas well as relatec 

0, fl H H 0, H Q 



and are omitted here. Note, 



however, that the width L of (both bright- and dark-types of) these localized excitations 
depends on the maximum amplitude po as L ■ 



'2\P/ {Q po)|; therefore, we retain that for a 
given amplitude, the (absolute value of the) coefficient ratio P/Q expresses the square width 
of the soliton, i.e. a pulse if PQ > and a hole if PQ < 0. Inversely, for a fixed width L, 
the quotient P/Q expresses the amplitude (height) of the solitary wave po- 

In figures El - El we have depicted the ratio P/Q as expressed by the relations (|38|). 



(jH)), expressed in units, say: Pq/Qq 



2u)'i 



the presence of negative dust (p < 1, see Fig. Et-), the soliton width is seen to bear 
lower values (with a maximum for higher K) as p decreases; therefore, an increase in the 
concentration of negative dust results in generally narrower excitations, but with a peak at 
higher wavenumbers K. Also, for a given K, the width is maximum for a certain value of p 
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(see Fig. Eb); the position of the maximum depends only shghtly on 9 but rather strongly 
on K (see Figs, lit-? b). Finally, for a fixed value of yU, the P/Q vs. K curve seems to have 
a maximum at = 7r/2; see Fig. El transverse modulation slightly favours higher soliton 
widths. This maximum moves to higher K with increasing dust (i.e. decreasing /i); cf. Figs. 

For positive dust (/i > 1), see Figs. IHl-IHl we have similar qualitative results, yet generally 
lower values. Once more, the angle variable does not seem to influence the soliton profile 
dramatically. 

VII. CONCLUSIONS 

In this paper, we have studied the 2D linear and nonlinear propagation of R-D-MHD 
waves in a uniform cold magnetoplasma composed of electrons, ions, and charged dust grains. 
The presence of immobile charged dust grains is responsible for the ion rotation and a new 
cutoff frequency (non-existing in an ordinary e-i plasma), which were reported by Rao in his 
classic paper j9|. The propagation of the modified dust magnetoacoustic waves is possible 
due to the finite ion inertia effect. The charged dust modifies the phase speed of the modified 
magnetosonic waves. Furthermore, we have considered the amplitude modulation of the R- 
D-MHD waves and have shown that self-interactions among waves result in the harmonic 
generation and the amplitude modulation of a carrier R-D-MHD wave. The wave envelope 
has been shown to be stable against perturbations in a wide range of physical parameter 
spaces. Finally, we have discussed the possibility of localized envelope excitations (mostly 
of the dark soliton type i.e. localized field dips propagating in the plasma) associated with 
the nonlinear R-D-MHD. 
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APPENDIX A: IST-ORDER PERTURBATION: DERIVATION OF THE DIS- 
PERSION RELATION AND IST-HARMONIC AMPLITUDES 



Consider the system: (fT^ - ^lU^ . which describes the evolution of S = 
E, B). By hnearizing around the equihbrium state So = 
[riefl, riifi, 0, 0, 0, Bo) viz. S = So + Si and assuming hnear perturbations of the 
form: Si = Si expi(kx — ut) = Si expi{kx + ky — ut), we obtain a new system of (hnear) 
equations for the perturbation amphtudes {Si)^: 

-iufie^i + ik^nefiUei) = 0, (Al) 

- iuni^i + ik{ni^oU-a) = , (A2) 

1 ^ 

Elx = — - Uel,y Bq , (A3) 

Ely = ^ Uei,x Bo , (A4) 
rui {-i uj )uii^x = ZiC (Eix + -Ui^yBo] = H — '-^-^ {ui^y - Ue,y) , (A5) 



c / c 



and 



(-. u )ua,y = e [Ely - h^^^B,) = (^,,,. - ^,,,) , (A6) 

Attc 

ikyBi = —— {ZiriifiUii^x - ?^e,oMei,x) , (A7) 
Attc 

ik^Bi = — {ZiriifiUii^y - nefiUei,y) , (A8) 



i kx Ei^y - i ky Ei^x = ^ (^ tu) , (A9) 



where only first harmonic terms were retained. Now, eliminating the electron velocity am- 
plitudes from ()A5|1 - ()A8|) (i.e. solving for Ugij in the latter two and substituting in the 
former), one immediately obtains: 

I uj Uiix + s\lc,i Uii^y = i\lc,i- Bikx 

Heft ^T^ene,Q 

iujUily — sVtci—^^—^ Uii X = i^ci-, Biky (AlO) 

(the ion cyclotron frequency Qc,i was defined in the text). Also, one may substitute from 
fTOjl . (jXH) into (jXni) in order to obtain: 

k ■ Uel = , (All) 

-Do 
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and, once more, use 



6|) to eliminate Uei in it: 



-Do 



k • Uii . 



(A12) 



Now, (jAlOb . b), ()A12|) form a closed system, with respect to Unj (j = x,y) and Bi. In 
specific, one may solve the latter for Bi and substitute into the former two; one thus obtains 
precisely the system of equations (j^ . along with the definitions mentioned in the text. 
On a more systematic basis, one may define the matrix: 

/ _n7, , n niu ^ _ oiu ^ ^ n n n n n \ 



L«(o;,k) 



iluo 







ilkyUefl 




















—iluj 








ilk^Uifl 


ilkyUifl 




















Bo 








c 














-Bo 














c 
















—iluj 


























—iloj 

















47re 

—ne,o 





4vre 7 
—^i^ifl 






















47re 
— T^e.O 

































tlJvy 




-i^ 



V 

(A13) 

which arises naturally by isolating the /— th harmonic terms (at every order u) in equations 
(fT^ - (j2ni)- For instance, for n = / = 1, the system on top of this Appendix is formally 



expressed as: Lq^"* sf"* = 0. Now, the condition DetL)^' = leads exactly to the dispersion 
relation (PHj) . while the solution of the system is given by ()26|) - (jH^ in the text. 



(1) 



APPENDIX B: 2ND-ORDER PERTURBATION: GROUP VELOCITY, 0-TH 
AND 2-ND HARMONIC AMPLITUDE CORRECTIONS 

For n = 2, / = 1, we obtain the system of equations Lq^'' S^^"* = Ri^"*, where 
Lq^'* was defined in flA13|) and R^^'' denotes the vector: {^c^i^^ — Ueflcll^\ Xc^2^^ — 
UiflCr^^\ 0, 0, Xc^i^\ Xcq^\ 0, —1, — Cq^^^)'^ dB[^^ /d(. The compatibility condition im- 
posed in order for a solution to exist, can be formulated as the constraint: Deth^'' = 0, where 
is the matrix obtained by substituting the m— th column in Lo^"* by Rp"*. Whichever 
the choice of m (= 1,2,..., 9), by solving the resulting equation, one readily obtains the 
definition of A as the group velocity Vg = du/dk^ (as defined in the text). One then obtains 
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the solution 5*1^- = c^f^^ dBi^^ / dC, for (8 of) the elements of 82^^"*, in terms of one of them, 



e.g. of S'j^g = Bi'\ Assuming, with no loss of generality, that BY' = 0, one obtains for the 
coefficients c^^^'' the expressions: 

cf^ = 



(2) 



cf' = 

4-) - - 



.(21) 



|z 5o e Zi uI q^u^ + nl i) - 2nefinifiZinl^ + nj^VLl ^Z'^ 



+ Hen ^Iri 



-icrrii riefl OJ + BQeZi^ZiUi Q — n^fi) sin 20 



2cmi nl^Q u [nefliu^"^ - Vlli) + rii^QVtl^Zi] 



-Bq e riin Z- 



4ne_o ^ ^c,Me,0 - ZiUifl) cos^ 9 



.(21) 



1 TT mi He TliS) <jJ (? k"^ 



„(21) 
^6 



^(21) 



.(21) 



1 



1 TT mf riefl rii o uj (? k'^ 



Bp 

'< n emf uIqUJ k'^ 



iBp Zj 
\ IT emf uIqUOC^ k"^ 



.2^ ^ _ ^^2 u2{^o^i\-^i^i'0 + ^7^ {nefi- ZiUi^oY] cos 29 

+cmi [—Bp ZiUi pck'^ + SnieriefiUj {n^fi — ZiUifi) sin 29^ 
i Bp Zi [rrii Uip k^ + Sir (rig o — ^i^i,o)^] sin 29 

+8ti ecrrii riefl oj {uefl — ZiUifi) cos 29^ 
I Uefi UJ k'^ + eZi —87iecminefiU!{nefi — ZiUifi) cos29 
—iBp Zi [niinifi k"^ + 87re^ (ng^o — ZiUi pY] sin 29 | 
crrii —Bp ck'^riifi Zi + 8TrienefiUj{nefi — ZiUifi) sin 2^ 

+Bp Zi [miriifi k'^ + 87re^ (rig o — ZiUifiY] cos2^^ | 

(Bl) 



For n = 2, / = 0, we obtain the system of equations Lg"-* Sq^"* = Rq^^ I-BI'^-'P (set / = 
in (imi) for 4°^), where R?^ is the vector (0, 0, -c^^\ c'i^\ ikyc'i^'>cf^^* + ^(4^^^ - 

r(^^)\ ^(11)^(11)* f^c,./(ll) m)N 4^e/7 J11)J11)* J11)J11)*A 4^e/7 Jll) Jll)* 

Ci C4 ), 0) + c.c. The 1st, 2nd and 9th equations are identically satisfied, so the cor- 
responding equations for n = 3, / = have to be "borrowed" . Combining them with the 
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remaining (3rd to 8th) equations here, we obtain 



n, 



(2) 
eO 



(2) 



«0 



'^e.O (20) |r(1)|2 
^9 

^ifi (20) |r(1)|2 
— C5 1^1 I , 



u 



(2) 

^^0 



(2) 



2u 



cose\B\ 



(1)|2 



C 



(20) 1^(1) 1 2 



(20) I d(1)|2 



c 



(20) I r(1)|2 
5/6 



(2) 



E = B 



(2) 



(B2) 



For n = 2, / = 2, we obtain a system of (9) equations in the matrix form: Lq S 



(2) e(2) 



R^'^ BY^ [set / = 2 in fXTa] for L^'O; the (lengthy) expression of the vector is omitted. 
Solving for the second-harmonic amplitudes Sj^^\ we obtain: 



n. 



(2) 



e2 



(2) 



^e,o ["•^i^»,oC^ + 47re^(ne,o - nj^pZ^)^] (1^2 _ (22) 



(2) 



e,a;2 



1 (22) ^(1)2 _ (22) ^(1)2 

c m'^ uj 



GnenlQu'^ + vr e q nj^o (— 4a;^ + Ql 



cos 9 



-I Hen iO 



Bq c Qc,i (3^6,0 + ZiUifl) + TT e [-2 o(2t^^ + fi^ •) 



2 n2 1 



sin 



J22) ^(1)2 



(2) 



3 3 

c u; 



e,?/2 



Svre^ riifi Zf B^{nefl - ni^oZif k 
jf^ci 67r e o cu^ + vr e q nj^o (-46^;^ + 

+ i ne,o Bock"^ Q.c,i (3ne,o + ^i^i^o) + tt e [-2 0(26^^^ + ^c. 



sin 6 



+ riefi riifl Zi^l i + Zl fi; 



2 n2 1 



COS 6 



- J22) o(l)2 

= C4 



19 



u 



(2) 



t,X2 



(2) 



E. 



(2) 



x2 



E 



(2) 



X 



37re3 nf -B^(ne,o - UifiZiY k 



cos^ 



+37r i e ne,o f^c,i (^^e.o - n-ifiZi) sm6\ = , 



X 



I So c A;^ -^i^i.o a; + TT e [n^ o(2a;^ + Q^,?) " 2ne,o Ui^o Zi^^^ + Zln^^fll^^ 



sin^ 



-Stt i e ne,o <^ (^^e.O - "-i.O-^i) COS 6' [> = cf^^ Sf ^ , 



c 



37re3 nj,o -B^(ne,o - ni^ZiY k 



X 



cos 9 



+67reng o<^^ sin 



m? c 



37re3 nj,o Zf Bl{n^fl - riiflZiY k 



X 



Bock'^ {SnloUJ^ + n^fsni^^^ZiVL\- - n\^^Zl^\^ + 

_ (22) ^(1)2 



sin^ 



e2 ni,o Z? -Bo ("-6,0 - ^if^Zif 



(B3) 
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Figure Captions 

Figure 1. 

The value of the coefficient Q is depicted against the dust parameter /x and the (normal- 
ized) wavenumber K. 

Figure 2. 

Soliton solutions of the NLS equation for PQ < (holes); these excitations are of: (a) dark 
type, (b) grey type. Notice that the amplitude never reaches zero in (b). These excitations 
represent electromagnetic field dips (voids) associated with the nonlinear R-D-MHD wave 
propagation. 

Figure 3. 

Negative dust; the (normalized) soliton width L (absolute value of P/Q) is depicted: (a) 
against wavenumber for ^ = and /i = 0.8, 0.7, 0.6, 0.5 (from top to bottom); (b) against 
the dust parameter /i, for K = 0.2 and 9 = 0°, 30°, 60°, 90° (from bottom to top). 

Figure 4. 

Negative dust; the (normalized) soliton width L (absolute value of P/Q) is depicted 
versus the dust parameter /x and the angle 6 for: (a) K = 0.2; (b) K = 0.5. 

Figure 5. 

Negative dust; the soliton width L (absolute value of P/Q) is depicted versus the 
wavenumber K and the angle 9 for: (a) fj, = 0.8; (b) fj, = 0.5. 

Figure 6. 

Similar to Fig. El for positive dust; the (normalized) soliton width (absolute value of 
P/Q) is depicted: (a) against wavenumber K, for 6 = and fi = 1.5,2.0,2.5 (from top 
to bottom); (b) against the dust parameter /x, for K = 0.2 and 9 = 0°, 30°, 60°, 90° (from 
bottom to top). 
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Figure 7. 

Similar to Fig. HI for positive dust; the (normalized) soliton width (absolute value of 
P/Q) is depicted versus the dust parameter /i and the angle 6 for: (a) K = 0.2; (a) K = 0.5. 

Figure 8. 

Similar to Fig. 13 for positive dust; the sohton width (absolute value of P/Q) is depicted 
versus the wavenumber K and the angle 9 for: (a) /i = 1.2; (b) /i = 1.5. 



22 




FIG. 1: 



23 





FIG. 2: 



24 





0.2 0.4 0.6 0.8 1 



FIG. 3: 



25 



Negative dust, K = 0.2 




Negative dust, K = 0.5 
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